Isotropic-nematic phase behavior of length polydisperse hard rods 
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The isotropic-nematic phase behavior of length polydisperse hard rods with arbitrary length 
distributions is calculated. Within a numerical treatment of the polydisperse Onsager model using 
the Gaussian trial function Ansatz we determine the onset of isotropic-nematic phase separation, 
coming from a dilute isotropic phase and a dense nematic phase. We focus on parent systems whose 
lengths can be described by either a Schulz or a 'fat-tailed' log-normal distribution with appropriate 
lower and upper cutoff lengths. In both cases, very strong fractionation effects are observed for 
parent polydispersities larger than roughly 50 %. In these regimes, the isotropic and nematic phases 
are completely dominated by respectively the shortest and the longest rods in the system. Moreover, 
for the log-normal case, we predict triphasic isotropic-nematic-nematic equilibria to occur above a 
certain threshold polydispersity. By investigating the properties of the coexisting phases across 
the coexistence region for a particular set of cutoff lengths we explicitly show that the region of 
stable triphasic equilibria does not extend up to very large parent polydispersities but closes off at 
a consolute point located not far above the threshold polydispersity. The experimental relevance of 
the phenomenon is discussed. 



I. INTRODUCTION 

Dispersions of highly anisometric rod- or platelike col- 
loidal particles exceeding a certain concentration are 
known to undergo a transition from an isotropic state 
(I), in which the particles are oriented in random di- 
rections to an aligned nematic state (N) J^, 2, 3J. The 
theoretical basis for this phenomenon has been estab- 
lished by Onsager in the 1940s. In a seminal paper 0, 
he showed that the phase transition can be explained 
on the basis of purely repulsive interactions between the 
particles. The basic mechanism behind the transition is 
a competition between orientional entropy which favors 
the isotropic state and the entropy effect associated with 
the orientation-dependent excluded volume of the aniso- 
metrical particles which favors the ordered nematic. 

A main characteristic of systems of (anisometric) col- 
loidal particles is their inherent polydispersity, i.e. the 
particles may differ in size and shape [Hi]. The issue 
of polydispersity and its effect on the interpretation of 
experimental results has already been addressed by On- 
sager in his original paper. Later on, extensions of the 
Onsager theory allowing for phase diagram calculations 
for bidisperse 0, H, 0, liol llll IT^ and tridisperse sys- 
tems [13 of hard rods as well as binary mixtures of hard 
platelets revealed a rich variety in behavior, most 
notably a widening of the coexistence region, a fraction- 
ation effect (i.e. segregation of the species among the 
coexisting phases), a reentrant phenomenon and, most 
interestingly, the possibility of a demixing of the nematic 
phase which may give rise to isotropic-nematic-nematic 
triphasic equilibria. 

So far, very few theoretical attempts have been made 



'Electronic address: ^ G.J.Vroege@chem.uu.nT| 



to study the isotropic-nematic phase behavior of truly 
polydisperse systems, i.e. systems characterized by a con- 
tinuous distribution in particle size, within the Onsager 
treatment. As to rodlike systems, there are some stud- 
ies in which the effect of polydispersity on the isotropic- 
nematic transition is accounted for using perturbation 
theories The major drawback of this approach, 

in which the effect of polydispersity is considered as 
a small perturbation to a monodisperse reference sys- 
tem, is its limited applicability; it is only justified for 
slightly polydisperse systems with very narrow length 
distributions. Although the perturbation approach qual- 
itatively predicts some generic features such as a broad- 
ening of the coexistence area and a fractionation effect, 
with the longest rod going preferentially into the ne- 
matic phase, other interesting phenomena which are ex- 
pected to occur at much higher polydispersities (in par- 
ticular polydispersity-induced demixing transitions) are 
not treated and therefore remain elusive. 

Solving the phase equilibrium conditions for systems 
with arbitrary size distributions is by no means trivial 
and requires considerable numerical effort. In particular, 
the presence of almost infinitely many components in a 
polydisperse system requires an equally large number of 
coexistence conditions to be solved simultaneously which 
obviously is a formidable task. Recently, a number of 
studies have appeared in which the fully polydisperse On- 
sager model, albeit in simplified form, was subjected to a 
numerical treatment. Speranza and SoUich p?7lll8j inves- 
tigated the model using the so-called 7'2-approximation 
which consists of truncating the angular dependence 
of the orientation-dependent excluded-volume term af- 
ter the first nontrivial term in an expansion in terms 
of even Legendre polynomials 'P2n- A remarkable out- 
come of these calculations is that for rod length distribu- 
tions with sufficiently fat tails (e.g. log-normal distribu- 
tions) triphasic isotropic-nematic-nematic equilibria are 
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predicted to occur in a small interval of polydispersities 
of the parent system. However, the simplified orientation 
distribution function (ODF) for the nematic phase per- 
taining to the '7^2-moder is only valid for the description 
of very weakly aligned nematic phases. The behavior 
predicted from this model should therefore be considered 
with some care, particularly in those regions where the 
fractionation effect is strong and the phase behavior is 
dominated by the effect of the longest rods in the sys- 
tem. The presence of very long rods in a nematic state 
may force the entire system into a strongly aligned ne- 
matic configuration, so that a more appropriate form for 
the nematic ODF is required in these cases. 

In this paper we use the Gaussian trial ODF approach 
to calculate the isotropic-nematic phase behavior of hard 
rod systems which can be described by either a Schulz 
or a log-normal length distributions with arbitrary poly- 
dispersities. The benefit of using the Gaussian Ansatz is 
twofold. First, all necessary integrals for the monodis- 
perse Onsager model are analytically tractable so that 
only numerical integrations over the length distributions 
need to be considered for the polydisperse case. Sec- 
ond, the Gaussian ODF allows for a qualitatively better 
description of highly ordered nematic states compared 
to the 7'2-approximation which makes it a suitable tool 
for describing polydisperse systems, particularly the ones 
with a 'fat-tailed' length distribution. While work on this 
subject was still in progress, Speranza and SoUich re- 
ported a numerical analysis of the exact Onsager model 
[l9|, i.e. using the numerically exact ODF. Also there, 
triphasic equilibria were predicted for both Schulz-type 
parent distributions and 'fat-tailed' log-normal forms. 
However, due to the numerical complexity of the prob- 
lem only the onset of nematic ordering from an isotropic 
reference phase was considered there so that no infor- 
mation could be obtained about the properties of the 
isotropic and nematic phases across the coexistence re- 
gion. Consequently, no conclusive insight could be gained 
as to whether the triphasic equilibria constitute a signif- 
icant part of the phase diagram. Within the Gaussian 
approach it is possible to access the coexistence region 
with only limited additional numerical effort. An impor- 
tant consequence is that it also enables us to study the 
triphasic demixing in more detail and to gain insight in 
the extent of its stability region. Although we cannot 
calculate the binodal curves for these equilibria (which 
locate the precise onset), we are able to localize the spin- 
odal points for the coexisting nematic phase across the 
two-phase region. The presence of these spinodal points 
indicates a local instability of the nematic phase which 
implies that a triphasic demixing transition must take 
place. 

This paper is structured as follows. In Sec. II we 
describe the Onsager theory generalized to polydisperse 
systems. The conditions for phase equilibria are outlined 
in general terms in Sec. III-A, and specified to the onset 
of isotropic-nematic phase separation and the full phase 
split situation in Sees. III-B and III-C respectively. In 



III-D some details of the truncated parent distributions 
we used in our calculations are given. The numerical 
results for the different parent distributions will be pre- 
sented and discussed in detail in Sees. IV and V. Finally, 
we summarize and conclude in Sec. VI. Technical as- 
pects of the subject are treated in several appendices. 
In App. A we provide details about the numerical pro- 
cedure we adopted to solve the coexistence conditions. 
In App. B we explicitly show that the Gaussian Ansatz 
yields the exact scaling results for parent distributions 
with infinite cutoff lengths and in Sec. C we establish a 
spinodal instability criterion for the nematic phase inside 
the coexistence region which is used to detect possible 
isotropic-nematic-nematic triphasic equilibria along the 
coexistence trajectory. 



II. POLYDISPERSE ONSAGER THEORY; 
STARTING EQUATIONS 

Let us consider a system of hard rodlike cylinders with 
equal diameters D but different lengths L , in a macro- 
scopic volume V. To characterize the rod lengths in our 
polydisperse system we introduce the relative rod length 
I — L/Lq (with Lq some reference rod length Lq) which 
is assumed to be a continuous variable. We may then 
take the limit Lq/D — > 00 (infinitely thin rods) at con- 
stant values for the relative lengths I. In the Onsager 
approach, the excess free energy describing the excluded 
volume interactions between the particles is truncated af- 
ter the second virial term. This approach can be shown 
to yield the exact free energy in the limit of infinitely thin 
rods. A generalization of the original Onsager model to 
include polydispersity leads to the following expression 
for the total Helmholtz free energy density / (in units 
ksT = /3-I) 



/ ^ 



bf3F 



c{l)[\nc{l) - l]dl+ / c[l)uj{l)dl 



c{l)c{l')ll' p{l,l')dldl' . 



(1) 



All irrelevant contributions linear in c arising from the 
standard chemical potentials of the particles are omitted 
since they only depend on the solvent chemical potential 
and the temperature. The concentrations c are rendered 
dimensionless by relating them to the orientationally av- 
eraged excluded volume per particle between two refer- 
ence rods, b = ttDLI/A, via c(l) = bN{l)/V where N{l)dl 
is the number of particles with relative length between I 
and I + dl. The density distribution over lengths c(Z) can 
be decomposed according to c{l) — cop{l), with p{l) the 
normalized length distribution (J dlp{l) = 1) and cq the 
total dimensionless rod concentration. 

The free energy Eq. (QJ consists of several entropic 
contributions. The first term represents the exact ideal 
free energy of the polydisperse system. The second term 
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contains parameter w as a measure for (the negative) of 
the oricntational entropy 



uj{l)= 'iJj{l,n)\ii[ATrij;{l,n)]dn, 



(2) 



where ip{0^ I) is the normaUzed ODF for species I describ- 
ing the distribution of the particles' sohd angle fi. In the 
isotropic state, all orientations are equally probable so 
that V'iso is simply a constant {l/47r) independent of I. In 
the nematic state, however, the ODFs are peaked func- 
tions (generally different for each species I), due to the 
fact that the rods are aligned along a nematic director. 
Note that u (and hence the orientation free energy) at- 
tains its minimum (uj = 0) in the isotropic state, whereas 
w > in the nematic state. 

The last term in Eq. describes the excess free 

energy which accounts for the particle interactions. In 
a second virial approximation, the interactions between 
hard particles may be expressed as an excluded volume 
entropy depending on the excluded volume between two 
particles. A measure for the average excluded- volume in- 
teraction between rods of relative length / and /' is given 
by the following angular average 



p{l,l') 



(3) 

Using the isotropic average ((|sin7(ri, f2')|)) = 7r/4 we 
obtain p{l,l') = 1 for the isotropic state. This indicates 
that the excluded volume (or packing) free energy is in- 
deed maximal in the isotropic phase but decreases as soon 
as the rods align to form a nematic phase. 

In this study we use Gaussian trial ODFs with varia- 
tional parameter a{l) to describe the angular distribution 
of rods with relative length / in the nematic state • The 
Gaussian Ansatz consists of supposing 



a(l) r 
a(l) r 



-ia(0(7r 



< 6i< 



Note that, 



< 6* < TT 

(4) 

due to the 



where a is now a function of / 
uniaxial symmetry of the nematic phase, the ODFs only 
depends upon the polar angle between the particle ori- 
entation vector and the nematic director. Inserting Eq. 
Q in Eq. ij^Jl and straightforward integration yields for 
the oricntational entropy 



uj{l) - lna(/) - 1. 



(5) 



For the excluded volume entropy in the nematic phase 
Pnem(^^') Only the leading order term of its asymptotic 
expansion for large a will be retained 



-3/2(0, a-3/2(?')' 




(6) 



Substituting Eqs. lO and © into Eq. (QJ and mini- 
mizing the free energy density with respect to the non- 
conserved oricntational degrees of freedom by means of 
a functional differentiation with respect to a{l) gives 



Sail) 



411 

a{l) 



1/2 



lc{l) 



l'c{l') 



2a2(?) 
1 1 



-1/2 



dV. (7) 



Applying the stationarity condition 5f /5a{l) = and 
some rearranging leads to the following self-consistency 
equation 



1 



a{l') 



-1/2 



dl' 



(8) 



Here, we have factorized the Gaussian variational pa- 
rameter function a(l) into a concentration-dependent 
part and a contribution a{l) only related, via Eq.®, to 
the normalized length distribution in the nematic phase 
p^^\l). Hence we write 



4^2 

a{l)=a{l)^. 

TT 



(9) 



showing that for all I the variational parameter a de- 
pends quadratically on cq just as in the monodisperse 
case |2y|. An approximate analytical solution to Eq. © 
valid for infinitely narrow distributions 

[3 (denoted by 

subscript S) can be obtained by substituting a delta func- 
tion p^^\l) = S{1 — 1) which gives 



"5(0 = ^(^8^2 + 1-1) 



(10) 



This result may be interpreted as a measure for the ne- 
matic alignment of a single rod with relative length I 
added to a nematic bulk system of monodisperse rods 
with reference length Lq. Eq. (|10|l shows that asil) and 
hence the order parameter ^20] P{1) 1 — 3/a{l) are in 
general, as we might have anticipated, increasing func- 
tions of the relative rod length, i.e. as{l) oc I for large 
Moreover, as(0) — which means that there is no order- 
ing for rods of zero length, as formally must be the case. 
However, it should be pointed out that rods with lengths 
close to zero must be excluded from our model because 
the normalization factors for the Gaussian ODFs in Eq. 
(@J do not allow for a correct description of isotropically 
distributed or weakly aligned species in the nematic state. 
For consistency reasons we must therefore introduce a 
lower limit (/min > 0) in all length distributions. 



III. I-N PHASE COEXISTENCE 

A. Equilibrium conditions for polydisperse systems 

The conditions for phase equilibrium are that the co- 
existing isotropic and nematic phases must have equal 
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chemical potential fj,{l) for all relative rod lengths I, as 
well as equal osmotic pressure 11. The chemical potential 
can be derived by functional differentiation of the free 
energy with respect to the length distribution c{l) 



/3/i(0 



5c{l)' 



(11) 



Using Eqs. (0) and © together with the isotropic values 
a; = and p = 1 we obtain 



/3Mi.o(0 = lnc'-'\l) + 2lc[''> 
/3/incm(0 = lncW(0+ln 



- 1 



(12) 



where ci denotes the first moment density following from 
the definition 



^ conik = J dll''c{l), k = 0,1,2, 



(13) 



Here, denotes the fc-th moment of the (normalized) 
distribution. The excess chemical potential for the nc- 
matic phase is given by 



/iW(/) = 23/2/ / dl'p^''\l')l' 



1 



1/2 



(14) 



and is independent of the concentration of the nematic 
phase. Similarly to Eq. © we can straightforwardly 

obtain an analytical solution for fiix \l) valid for near 
monodisperse distributions by substituting p^^)(/) = 
6 {I — 1), which yields the following scaling result 



1 



as{l) 



+ 1. 



(15) 



Using Eq. ifTUI) it follows that oc for very large 

I. This means that the excess chemical potential (i.e. the 
reversibel work required to insert a single rod of length 
I in a nematic system of reference rods ) increases with 
length, which is consistent with intuition. The osmotic 
pressure can be written in terms of the chemical potential 
and the free energy via 



(16) 



bpU=-f + P j dlc{l)^iil), 
which immediately yields for the isotropic phase 



b(3U 



ISO '-'0 



(17) 



For the nematic phase Eq. Ijltil) this formally gives 

6/3n„em - 4"^^ + /if ^ (18) 

However, this result can be simplified considerably by 
noting that the volume fraction of the average excluded 



volume (per particle) in the nematic phase is a constant, 
namely 



N 



V 



c[''H{U'p{l,l'))),,^2. (19) 



The brackets denote averages over the normalized length 
distribution. This result, which is due to Odijk gen- 
erally holds for both monodisperse and polydisperse sys- 
tems, independent of their composition. From the free 
energy Eq. it then follows that /, 



2cq ' so that 



the osmotic pressure of the nematic phase reduces to 



3c! 



(TV) 



(20) 



like for a monodisperse system We can now state 

the conditions for the coexistence between the isotropic 
and nematic daughter phases into which a parent phase 
(henceforth denoted with superscript 0) with length dis- 
tribution c^-^\l) is assumed to have split. From Eq. H12|l . 
equality of chemical potentials of both phases is obeyed 
exactly if the distributions in the phases have the follow- 
ing form 



a = I,N (21) 



where W{1) = exp[/3/i(Z)] must be a function common to 
both phases, since iJ,'^^\l) = /i*-^-'(Z) = A*(0- The func- 
tions ^^"-'(0 are given by 

^^'Hl) = -2lc['^ 

CW(0 = (l-lnl)-21ncr-ln5(0-A(^)(0. 

(22) 

Furthermore, conservation of matter requires 

c^°\l) = {l-j)c^''>{l)+jc^^'>{l), (23) 

where 7 denotes the fraction of the system volume oc- 
cupied by the nematic phase. Using Eq. H23|l . we can 
express W{1) in terms of the parent distribution c'^^^l) 
which gives 



(0 



.(0) 



(0 



exp 



(l-7)expKW(0]+7expKW(/)]- 



(24) 



These functions represent the equilibrium rod length dis- 
tributions for the coexisting phases. The phase equilibria 
can now, in principle, be obtained by solving a set of self- 
consistency equations for the moment densities of both 
phases and for the functions a{l) and p,ix\l) pertaining 
to the nematic phase. These equations will be worked 
out below for a specific situation, namely at the onset of 
isotropic-nematic phase separation. 
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B. The onset of I-N phase separation; cloud and 
shadow curves 



In this section we aim at locating the onset of isotropic- 
nematic phase separation indicated by so-caUed cloud 
and shadow points. A cloud point marks the density 
where a parent phase starts to split off an infinitesimal 
amount of a new coexisting phase, called the shadow 
phase. Accordingly, at the isotropic cloud point only an 
infinitesimal amount of nematic phase (shadow phase) 
has emerged and so the distribution of the isotropic 
phase is only negligibly perturbed away from the parent. 
Hence, for the isotropic cloud point we may set 7 = in 
Eq. ^ so that. 



(0 



.(0) 



(0, 



(25) 



implying that the distribution in the isotropic phase at 
the cloud point is equal to the parent distribution and 



^0 



.(0) 



ug . The (normalized) rod distribution 

in the nematic shadow phase (with density Cq^'') is now 
given by 



(0 



„(^) 



exp 



(0 



N — 



■ exp 



2c(");-A(f(0l, (26) 



where /Cjv = '^^c^o'^ ^'^'^ p^'^^l) the normalized 

parent distribution. Note that Eq. H2()|l is an implicit 
expression for p^^\l) because it still depends on the un- 
known functions for the variational parameter a{l) and 

the excess chemical potential for each species in 

the nematic shadow phase. Explicit solutions for these 
functions can be obtained by substituting Eq. (|26|l into 
Eqs. JHl and 1)14(1 and numerically solving the resulting 
self-consistency equations. 

The concentrations of the isotropic cloud phase and 
the coexisting nematic shadow are found by imposing 
the normalization condition for the distribution in the 
nematic shadow phase. 



pW(0 = i, 



and the condition of equal osmotic pressure 



(27) 



(0) 



(0)x2 



+ (cr) 



(28) 



Using this simple equation to eliminate e.g. 
may conveniently combine Eqs. 1)27(1 and ((28(1 into one 
self-consistency equation for the concentration of the 
isotropic cloud point, which we can solve in an iterative 
fashion. However, since a{l) and fiii^\l) also depend on 
Co°^ [via p^^^l), Eq. this equation has to be solved 

along with the coupled set of self-consistency equations, 
Eqs. ijHl) and l(14() . so that we end up with a set of three 



coupled nonlinear equations. Obviously, solving this set 
is not a trivial task but requires some numerical effort. 
For this reason we have devoted an Appendix A to this 
issue in which we describe some details of the numerical 
procedures adopted in this study. 

We can now perform a similar analysis to obtain ex- 
pressions for the nematic cloud point and the associ- 
ated isotropic shadow point, which locate the onset of 
I-N equilibrium coming from a dense nematic parent 
phase. Since the latter now coexists with an infinitesimal 
amount of an isotropic shadow phase we may set 7 = 1 in 



(N) 



Co"-*. The (nor- 



^ so that c(^'(0 = c(°)(/) and c}^ 
malized) rod distribution in the isotropic shadow phase 

(with density Cq^-*) is then given by 



.(0) 



(0 



exp 



ICip^°\l)&{l)exp \~2c^pl + ili^\l)\ ,(29) 



where ICj = A{c'jj '')'^ /necQ^K Since the functions a{l) and 
ij'C7i\l) correspond to the nematic parent phase we may 
substitute p(^)(0 = P^°Hl) in Eqs. © and Con- 
sequently, as the normalized parent distributions p*^°^(^) 
have a predescribed form and do not depend on any con- 
centration, the coupled self-consistency equations need 
to be solved only once for a given p'"' (Z) . The forms of 
the parent distributions will be specified in Sec. III-D. 

Once the solutions for d{l) and fli^\l) have been 
obtained, the concentrations of the nematic cloud and 
shadow point can be calculated by requiring self- 
consistency for the zeroth moment (normalization con- 
dition) and the first moment of the isotropic shadow dis- 
tribution, i.e. 



P^'Hi) 



= 1, 



(30) 



These conditions can, similarly to the previous case, eas- 
ily be rewritten into a set of consistency relations for the 
densities of the nematic cloud and shadow points, which 
we can solve iteratively (see also Appendix A). 



C. Inside the coexistence region 

We will now focus on the coexistence region between 
the isotropic and nematic cloud points, where both 
phases coexist in noninfinitesimal amounts, i.e. < 7 < 
1. According to Eq. ((24(1 . the equilibrium length distri- 
butions in the coexisting phases are then given by 



(0 



(0 



{AC[ci(0,A^x^(0]} + (l-7)' 



7 exp • 
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(0 



(1 



7)exp{-Ae[ci(/),Mi^'(/)]}+7 



•(31) 



with A^(0 = ^(^HO - ^^^HO, given by Eq. 
Note that both distributions are now different from the 
parental one. Solving the coexistence problem is done 
in a similar way to the one described in Sec. III-B for 
the isotropic cloud and shadow points. From an exper- 
imental point of view, we are only interested in results 
located on so-called physical dilution lines along which 
the overall system number density (c^'^) is changed (by 
e.g. adding or evaporating solvent) while the overall com- 
position of the species (p*^"^ (l)) remains fixed. The parent 
distributions will be specified below. 




D. Parent distributions 

The numerical method described in the previous sec- 
tions allows us to calculate the isotropic-nematic phase 
diagram for in principle arbitrary parent distributions. 
In our study we specify two types of distributions. The 
first one is the Schulz distribution which has the form 



= Nl'exp[-{z + l)l], 



(32) 



with normalization factor N. In order to exclude rods 
with potentially zero length we truncate the distribu- 
tion at some lower cutoff length Imin- Henceforth we fix 
^min = 0.01. For calculational purposes (see Appendix A) 
we must also have some finite cutoff length Imax at large 
I. Of course, introducing finite cutoff lengths is also rea- 
sonable from a physical standpoint. The first and second 
moment (defined as = {I''), k = 1,2) of the Schulz 
distribution are mi — 1 and m2 = (z + 2)/(z + 1) only 
for the unbounded case. However for finite cutoff lengths 
the moments will deviate from these values. Although 
the corrections are generally small, in particular for large 
^max, they cannot be neglected. Therefore, we choose to 
calculate all relevant moments of the parent distribution 
numerically via — ji""""^ l^p^^^dl. The polydispersity 
a is defined as 



{If 



1, 



and would yield = (1 + z)^^ for the unbounded Schulz 
distribution. 

The second distribution we consider is the log-normal 
one. The 'fat-tailed' log-normal distribution decays much 
slower at large I than the Schulz one and therefore pos- 
sesses a significantly larger contribution of long rods. The 
log-normal distribution reads 



exp 



In/ 



2w^ 



(34) 



For the unbounded log-normal distribution, w is directly 



FIG. 1: Concentrations of the isotropic and nematic cloud 
phases (solid lines) and the corresponding shadow phases 
(dotted lines) plotted against (on the vertical axis) the poly- 
dispersity (T of a Schulz parent with Imin = 0.01 and Zmax = 
100. The isotropic cloud curve is the one with the lowest con- 
centration. In the monodisperse limit {a = 0) the isotropic 
cloud point meets the shadow of the nematic cloud point and 
vice versa. The thin solid lines are the limiting curves for 
'max — > oo, given by Eq. H42|l in Appendix B. 



the parameter fx is chosen such that mi = 1, giving fi = 
—w/2. The second moment is then given by m2 = 1 + 
Also here, truncation of the distribution at some finite 
values /mill and /max leads to deviations for which we 
correct numerically. 



related to the polydispersity via w 



ln(l 



IV. RESULTS FOR THE ONSET OF I-N PHASE 
SEPARATION 

A. Schulz distributions 

In Fig. 1 to 5 we have depicted the results for a Schulz 
parent distribution with cutoff lengths /min — 0.01 and 
^max = 100. The curves describing the densities of the 
cloud and shadow phases are shown in Figs. 1 and 2. 
(33) A striking broadening of the coexistence gap can be 

detected, mainly due to a dramatic increase of the con- 
centration of the nematic cloud phase. In Fig. 2 we 
see that the volume fraction of the nematic cloud phase 
increases by several orders of magnitude at cr > 0.4. Al- 
though the nematic shadow curve crosses the correspond- 
ing isotropic cloud curve at a ~ 0.5 in Fig. 1, the vol- 
ume fraction (and hence the mass density) of the nematic 
shadow remains higher than that of the isotropic cloud 
phase throughout the phase diagram as we see in Fig. 2. 
Fig. 3 shows the extent of fractionation, i.e. the reparti- 
tioning of the long and short rods, among the coexisting 
phases at the onset of phase separation. A marked fea- 
ture is the rapid increase of the average length in the ne- 
^) and matic shadow around some 'transitional' polydispersity 
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FIG. 2; Similar to FIG. 1. but with the logarithm of the 
scaled volume fraction 4)Lq / D plotted versus the parent poly- 
dispersity a on the vertical axis. Note the dramatic increase 
of the volume fraction of the nematic cloud phase. 
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FIG. 3: Average length (Z) = ml in the isotropic and nematic 
shadow phases relative to the average length mi in the cloud 
phase plotted against (on the vertical axis) the polydispersity 
cr of a Schulz parent with Zmin = 0.01 and Zmax = 100. The 
inset shows the relative average length in the isotropic shadow 
phase (corresponding to the nematic cloud point). 



at — 0.5. This indicates that the nematic phase becomes 
preferably populated by the longest rods in the system. 
Note that there is a similar effect in the isotropic shadow 
phase around at — 0.7 where the shortest rods completely 
dominate the isotropic shadow phase at higher polydis- 
persities. The same effects are reflected somewhat clearer 
in Fig. 4 showing the evolution of the polydispersity of 
the shadow phases. At a = at, the polydispersity of 
the shadow phases show a kink. The strong decrease at 
higher a is due to the effect that the shadow phases be- 
come more and more enriched in either the longest or the 
shortest rods in the distribution. The dramatic change 
of the composition across the cr-range is shown explicitly 



FIG. 4: Similar to FIG. 3. but with the polydispersity of 
the isotropic and nematic shadow phases as (relative to the 
parental one a) plotted against a. Note the kinks in the 
isotropic and nematic branches around a — 0.65 and 0.5, 
respectively. 
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FIG. 5: Normalized length distributions in the nematic 
shadow phases p'^^ (I) for various parent polydispersities a 
of the same Schulz parent as in the previous figures. 



in Fig. 5 where we have depicted the normalized length 
distributions in the nematic shadow phase for various 
a. A similar picture is obtained for the distributions in 
the isotropic shadow phase (not shown here) but with 
the peak of the distribution shifting rapidly towards the 
lower cutoff length Zmin = 0.01. 

In summary, we can state that there are two fraction- 
ation regimes for the onset of phase separation. First, at 
low polydispersities {a < at) moderate fractionation is 
observed and the shadow phases are mainly populated by 
rods with slightly higher (or lower) than average length. 
Second, at higher polydispersity {a > at) strong fraction- 
ation occurs such that the shadow phases are completely 
dominated by the longest (or shortest) rods in the distri- 
bution. In a small interval around a = at the location of 
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FIG. 6: Concentrations of the isotropic and nematic cloud 
phases (solid lines) and the corresponding shadow phases 
(dotted lines) plotted against (on the vertical axis) the poly- 
dispersity a of a log-normal parent distribution with the same 
lower cutoff length Imin = 0.01 but two different higher cut- 
off lengths, imax ~ 10 and 100. The isotropic cloud curve is 
the one with the lowest concentration. The nematic cloud 
curve and the associated isotropic shadow curve are insensi- 
tive to the value of imax on this scale and therefore the results 
for iniax ~ 10 and 100 overlap. At the kink in the isotropic 
cloud curve the isotropic cloud phase coexists with two ne- 
matic phases differing in composition. 



the peak of the length distribution shifts rapidly, upon 
increasing a, from a value slightly different from one (per- 
taining to the low-cr regime) to a value close to the the 
cutoff length (corresponding to the high-cr regime) . 

So far, we have only discussed the results for a single 
set of cutoff lengths. Although the results for different 
cutoff lengths can be significantly different, particularly 
in the 'cutoff-dependent' regime ^3 c > cTf, the global 
phase behavior remains qualitatively the same. There- 
fore we conclude that the aforementioned fractionation 
scenario generally holds for any Schulz parent with suffi- 
ciently extreme cutoff lengths ( Imin *C 1 and Imax ^ !)• 
An interesting limiting case however is the behavior for 
Schulz parents with infinitely large cutoff lengths. For 
this specific case, we could obtain simple scaling rela- 
tions which describe the global behavior of the nematic 
shadow in the limit of an unbounded Schulz parent, i.e. 
^max — *■ oo. The scaling analysis, worked out briefly in 
Appendix B, is closely related to a more elaborate anal- 
ysis presented in Ref. [l^ for the exact Onsager model. 
In particular, we show in Appendix B that the Gaussian 
Ansatz must yield the exact high-cutoff scaling relations. 
The reason for this is that our high-cutoff scaling form 
for p(^)(Z) (for the nematic shadow phase) is completely 
analogous to the exact scaling form. 

Finally, we remark that we do not observe a real jump 
in the shadow curves (and a kink in the associated cloud 
curves), as found in Ref. The presence of a jump 

in the shadow curve indicates that, at some point, the 




"log((^L„/D) 

FIG. 7: Similar to FIG. 6. but with the logarithm of the 
scaled volume fraction (jiLo / D plotted versus the parent poly- 
dispersity a on the vertical axis. The main graph shows the 
results for the isotropic cloud and the nematic shadow phase. 
Note that the jump in the shadow curve corresponds to a kink 
in the associated cloud curve. The inset shows the results for 
the nematic cloud and its isotropic shadow phase. 



isotropic cloud phase coexists with two different nematic 
shadow phases and that a region of stable triphasic equi- 
libria developes. Therefore we conclude that, within the 
Gaussian Ansatz, the Schulz form does not give rise to a 
three-phase separation, at least up to /max = 100. 



B. Log-normal distributions 

The results for the log-normal case are presented in 
Fig. 6 to 8. The cloud and shadow curves shown in 
these figures correspond to log-normal parent distribu- 
tions with the same lower cutoff length /,nin = 0.01 but 
two different higher cutoffs /max — 10 and 100. We see 
that the phase behavior is globally the same as for the 
Schulz case; there is a generic broadening of the bipha- 
sic region (Fig. 7) and a very pronounced fractionation 
effect, particularly for the isotropic cloud and nematic 
shadow phases, as visible in Fig. 8. At low polydispersi- 
ties the distribution in the nematic shadow phase is very 
similar to the parent distribution (albeit with a higher 
average length). At higher polydispersities, however, we 
enter a regime characterized by extreme fractionation, i.e. 
the nematic shadow phase is dominated by the longest 
rods in the distribution. Remarkably, we do not see a 
similar transition in the isotropic shadow phase in Fig. 
8, as we did in the Schulz case. This implies that the frac- 
tionation from the nematic cloud phase is much weaker 
for log-normal distributions than for Schulz ones. 

A crucial difference with the previous results is that 
the transition between the regimes occurs discontinu- 
ously, that is, at the transition polydispersities at the 
isotropic cloud curves show a kink and the associated 
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FIG. 8: Average length mf in the shadow phases relative to 
the parental one mi for log-normal parent distributions with 
cutoff lengths Zmax ~ 10 and 100 as a function of the parent 
polydispersity a. The inset shows the relative average length 
in the isotropic shadow phase (corresponding to the nematic 
cloud point). Also here, the results for Zmax ~ 10 and 100 
overlap. 



nematic shadow curves exhibit a jump. Precisely at the 
kink, the isotropic cloud phase coexists with two different 
nematic phases, one containing mostly rods with slightly 
higher than average length (denoted by Ni) while the 
second one (Njj) is dominated by the longest rods in the 
distribution. Therefore, this special point marks the be- 
ginning of a stability region for I — Nj — Njj triphasic 
equilibria for log-normal distributions. In Fig. 8 we see 
that the position of the kink (in terms of at) rapidly shifts 
to lower polydispersities as ^max increases. From this, we 
anticipate that the triphasic equilibrium sets in at almost 
zero polydispersity (i.e. near monodisperse systems) for 
very large cutoff lengths. 

Like for the Schulz case, we can obtain information 
about the global phase behavior for parent distributions 
at infinite cutoff lengths /max from the high-cutoff scal- 
ing results, shown in Appendix B. The most important 
outcome is that the concentrations of the isotropic cloud 
and nematic shadow phases go to zero for large cutoff 
length rather than approaching asymptotic forms such 
as in the Schulz case. Furthermore, it is shown explicitly 
that the fractionation between the isotropic cloud and 
nematic shadow phases is stronger than for the Schulz 
case. 

So far, we have only looked at the onset of phase equi- 
librium by analyzing the properties of the cloud and 
shadow phases. The next step is to explore the coex- 
istence region in more detail. An intriguing issue is to 
verify the region of stability for the isotropic-nematic- 
nematic triphasic region for the log-normal case. This 
will be dealt with in the next section, where we discuss 
the phase diagram for a log-normal parent with cutoff 
length Zniax = 10 in more detail. 



m. 
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FIG. 9: Average rod length mi in the coexisting daughter 
phases (solid lines) plotted versus the concentration Cq"' of the 
parent across the coexistence region for a log-normal parent 
with Imin = 0.01 and /max = 10 at fixed polydispersity a — 
0.4. The curve for which mf > 1 is the nematic branch. The 
dotted line corresponds to mi — 1 for the parent phase. The 
inset shows the (logarithm of the) fraction 7 of the system 
volume occupied by the nematic phase. Note that the amount 
of nematic phase is extremely small in the region Cq"' < 3. 




FIG. 10: Evolution of the scaled volume fraction (pLo/D in 
the coexisting daughter phases across the coexistence region 
for the same parent. The nematic branch is the one with the 
highest volume fraction. The dilution line is indicated by the 
dotted line. 



V. INSIDE THE I-N COEXISTENCE REGION 

Across the coexistence region the equilibium length 
distributions of the coexisting phases, which originate 
from a parent phase with a prescribed distribution 
p^"^(/)), change continuously as the overall density of the 
parent is Cq"' is changed. In the actual calculations how- 
ever it is more convenient to impose the fraction 7 occu- 
pied by the nematic phase rather than Cq""* and calculate 
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FIG. 11: (a, left) Relative polydispersity af/cr of the coexisting daughter phases across the coexistence region for the same 
parent. By definition, the dilution line (dotted) is given by of /a = 1. (b, right) Plot of the normalized length distributions of 
the nematic phase across the coexistence region corresponding to the positions A through E in (a). 



the corresponding densities self-consistently. In Fig. 9 to 
12 we show the properties of the coexisting isotropic and 
nematic phases for a log-normal parent distribution with 
polydispersity a = 0.4. Fig. 9 shows that the average 
length in the nematic daughter phase decreases rapidly 
in the regime Cg*^^ < 3 whereas only weak changes are 
notable at higher Cq"'' . Furthermore, we see that the vol- 
ume occupied by the nematic phase is extremely small in 
this regime. The same feature is observed in the volume 
fraction representation in Fig. 10. In particular, the non- 
monotonicity of the nematic branch is reflected somewhat 
clearer here. The fact that the isotropic branch runs ex- 
tremely close to the dilution line for Cg°'' < 3 indicates 
that the fraction of nematic phase formed must indeed 
be very small. The rather exotic oscillations in the be- 
havior of the polydispersities of the daughter phases in 
Fig. 11(a) reflect the dramatic change of the shape of the 
length distribution in the nematic phase in the first part 
(7 < 10^^) of the dilution trajectory as shown in Fig. 
11(b). Note that the distribution of the nematic shadow 
phase is in fact bimodal, with a small peak around / = 1 
and a much larger one at I = ^max- When the overall den- 
sity is increased the second peak shifts to lower values of 
I and eventually coincides with the first peak. When the 
overall density has reached about Cg" — 4 (corresponding 
to a nematic phase volume fraction of about 10 %) the 
distribution of the nematic phase resembles the parental 
one, albeit with a slightly higher average length. 

In Fig. 12 we have plotted the evolution of the aver- 
age length for a parent with polydispersity a = 0.3. A 
peculiar behavior is observed, which is also reflected in 
Fig. 13 where the coexistence pressure is plotted versus 
7. Clearly, there is a region where the pressure decreases 
as function of 7 which suggests an instability (or van 
der Waals) loop indicating a possible destabilization of 
the nematic phase. In Appendix C we explicitly show 
that the local extrema in the osmotic pressure in Fig. 13 
represent spinodal points for the nematic phase which in- 



dicate that the coexisting nematic phase indeed becomes 
locally (and hence globally) unstable. In the region be- 
tween the spinodal points (where 5{bpiV)/5^ is negative), 
the coexistence between the isotropic and a single ne- 
matic phase also becomes unstable such that a triphasic 
isotropic- nematic- nematic (/ — Nj — Nn) demixing oc- 
curs. 

It should be stressed that the actual onset of the three- 
phase separation is marked by binodal points which we 
have not located in this study, except for the kink at at- 
In general, binodal points are located at a lower concen- 
tration than the spinodal points so that the demixing 
usually occurs well before the point where the system 
becomes locally unstable. This becomes clear in Fig. 
14 where we show the details of the phase diagram in 
the vicinity of the kink including the spinodal curves (in 

terms of (cq"'')) for the coexisting nematic phase. At the 
kink at — 0.264 the three-phase separation sets in right 
at the isotropic cloud point but the spinodal points are 
located at higher concentrations. An important feature 
in Fig. 14 is the presence of a high-cr consolute (or criti- 
cal) point at (T = 0.373 ±0.001 where the spinodal curves 
meet. This means that the region of stable triphasic equi- 
libria does not extend up to large parent polydispersities 
but closes off at the consolute point which also consti- 
tutes an endpoint of the binodal curves corresponding to 
the triphasic equilibria. 

Consequently, we can distinguish three regimes in the 
phase diagram depicted in Fig. 14. First, below the 
kink there is a common isotropic-nematic phase separa- 
tion involving a moderately fractionated A'j-phase. Sec- 
ond, in the region between the kink and the consolute 
point 0.264 < a < 0.373 a strongly fractionated Njj- 
phase (containing the longest rods) splits off initially at 
the isotropic cloud point. At higher concentrations (when 
/ and Nji coexist in finite amounts) a second nematic 
A'j-phase is formed and an / — Nj — Nu triphasic equi- 
libium develops. Upon slightly further concentrating the 
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FIG. 12; Average rod length mf in the nematic phase across 
the coexistence region for a log-normal parent with the same 
cutoff lengths at cr = 0.3. The inset shows a hysteresis loop 
indicating that the nematic phase becomes locally unstable. 




FIG. 13: Coexistence pressure across the coexistence region 
for the same parent. In the region where 5(6/311) /57 < the 
nematic phase is locally unstable with respect to a nematic- 
nematic demixing. 

sample the Nij-phase disappears and a regular / — Ni 
biphasic equilibrium is recovered. The third region is lo- 
cated at parent polydispersities above the consolutc point 
a > 0.373. Upon concentrating the isotropic phase a 
strongly fractionated nematic phase (reminiscent of the 
Nil phase) is formed initially but the composition of this 
phase evolves gradually towards a A^/-type nematic phase 
as the biphasic region is crossed. 

VI. DISCUSSION AND CONCLUSIONS 

We have numerically investigated Onsager's second 
virial theory for polydisperse hard rods within the Gaus- 
sian Ansatz. The onset of isotropic-nematic phase sep- 
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FIG. 14: Phase diagram for log-normal parent distribution 
with cutoff lengths Zmin = 0.01 and Zmax = 10. The thick 
curves delimit the spinodal instability region for the coexisting 
nematic phase (in terms of the parental concentration Cq"') 
indicating a triphasic I — Nj — Njj demixing. 



aration is obtained from the cloud and shadow curves, 
which delimit the coexistence region. Within the same 
numerical framework, we could also explore the proper- 
ties of the coexisting phases across the coexistence re- 
gion. In this paper, we focussed on systems of poly- 
disperse hard rods whose lengths can be described by a 
Schulz or a log-normal distribution. The basic difference 
between these two forms is that the fat-tailed log-normal 
one contains a significantly higher fraction of longer rods. 
For numerical and consistency purposes we truncated the 
distributions at both ends at sufficiently low and high 
cutoff lengths. Using truncated distributions is also jus- 
tifiable from an experimental standpoint. For parent dis- 
tributions of the Schulz type the phase diagram contains 
two fractionation regimes. First, at low parent polydis- 
persities moderate fractionation occurs and the average 
rod lengths in the isotropic and nematic shadow phases 
are not much different from the average length in the 
parental cloud phase. Second, at higher parental a the 
fractionation effect is extremely strong and the isotropic 
and nematic shadow phases are completely dominated by 
respectively the longest and the shortest rods present in 
the system. 

For the exact Onsager model, Speranza and SoUich 
[l9l | very recently predicted a kink in the isotropic cloud 
curve (and a jump in the corresponding nematic shadow 
curve) for Schulz parents with cutoff lengths Zmax > 50. 
The presence of this kink indicates a region of stable 
isotropic-nematic-ncmatic triphasic equilibria. Here we 
do not find any indication for such a three-phase separa- 
tion for Schulz parents at least up to Zmax = 100. The 
discrepancy may be due to the Gaussian Ansatz, which 
implies that the ODFs are not represented by their cor- 
rect equilibrium forms. Moreover, the Schulz form might 
be a borderline case since its tail is too modest to induce 
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a strong demixing but too 'fat' to suppress it completely 
so that the presence of a kink in the isotropic cloud curve 
depends quite sensitively on the precise representation of 
the ODF. 

Although the Gaussian ODF is not a solution of the 
exact stationarity condition for the ODF, it does satisfy 
the exact high-density scaling relation |23|. This means 
that the properties of highly ordered nematic states are 
described very well by the Gaussian form. In fact, the 
description becomes exact for infinitely aligned states. A 
manifestation of this is the osmotic pressure for the ne- 
matic phase, Eq. I|20|l . which is the exact high-density 
result [231 ■ Consequently, for our polydisperse systems, 
we expect the Gaussian Ansatz to work increasingly well 
both for highly concentrated nematic phases and nemat- 
ics that are dominated by the longest rods. In both cases, 
the nematic alignment of all species is expected to be very 
pronounced such that the use of the scaling ODF (for all 

1) is justified. To verify this notion, we have plotted the 
variational parameter as a function of length, for both 
the nematic shadow phase and the nematic cloud phase 
corresponding to a Schulz parent in Fig. 15. Since the 
Gaussian ODF is expected to be the least correct for the 
shortest rods (which show the weakest alignment) , we fo- 
cus on the interval Zmax < 1. In order for the results to be 
self-consistent, the alignment must be strong enough and 
hence the variational parameter must be sufficiently large 
(say ^°loga > 1) for all lengths. Fig. 15 (a) shows that 
this is not entirely the case; in the regime of low fraction- 
ation the shortest rods (with / < 0.4) are not sufficiently 
aligned by the longer rods so that the Gaussian descrip- 
tion fails here. In the regime cr > 0.5, which is physically 
the most interesting one, the ordering of the short rods is 
much higher due to the presence of very long rods in the 
nematic shadow phase, and the Gaussian Ansatz is fully 
justified. Similarly for the nematic cloud phase in Fig. 
15 (b) we see that the shortest rods are not well repre- 
sented by the Gaussian ODF at low a but much better at 
a > at where the variational parameter increases several 
orders of magnitude due to a dramatic increase of the 
concentration of the nematic cloud phase (see also Fig. 

2) . Obviously, for any rod length significantly larger than 
^min the Gaussian ODF works very well because a gen- 
erally becomes extremely large for any a. Therefore we 
conclude that given the fact that the composition of the 
nematic phases is dominated by the longest rods, particu- 
larly in the physically relevant cutoff-dependent regime, 
the Gaussian description is an appropriate tool in our 
study. 

We now turn to the log-normal case. The fractiona- 
tion scenario we observe there is qualitatively the same 
as the one for the Schulz case; weak fractionation occurs 
at low a where the distributions in the shadow phases are 
reminiscent of the parental one, but dramatic segregation 
effects take place above some threshold a, particularly be- 
tween the isotropic cloud and the nematic shadow phase. 
A crucial difference with the results for the Schulz case 
however is that the transition between the two fraction- 




FIG. 15: Behavior of the Gaussian variational parameter a{l) 
for the shortest rods of a Schulz parent with /min ~ 0.01 and 
'max = 100 (see also Fig. 1). (a, top) Results for the nematic 
shadow phase at various a. (b, bottom) Same for the nematic 
cloud phase. For comparison the result for the monodisperse 
system (a — 33.4) is indicated by the dotted line. 

ation regions shows a discontinuity at some a = at- At 
this point, the isotropic cloud curve shows a kink which 
corresponds to a jump in the nematic shadow phase. The 
jump indicates that an isotropic cloud phase must co- 
exist with two different shadow phases, one containing 
mostly rods of slightly higher than average length (the 
iV/-phase) and the other one predominantly containing 
very long rods (the iV//-phase) . The kink also marks the 
beginning of a region of stable I — Nj — Njj triphasic 
equilibria. For log-normal distributions with a moderate 
cutoff- length l^nx = 10 we found that the triphasic region 
opens up at a fairly low polydispersity at = 0.264. This 
value will decrease for larger cutoff lengths and eventu- 
ally go to zero when Ima.^ approaches infinity. This indi- 
cates that adding a very small fraction of long rods to a 
weakly polydisperse system of much shorter rods can al- 
ready induce a three-phase demixing. A similar effect is 
observed in binary mixtures of long and short hard rods 
with sufficiently large length ratios ^3 ■ 

By numerically analyzing the spinodal instability cri- 
teria for the nematic phase across the coexistence re- 
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gion for the case Zmax = 10 we have shown that the 
triphasic area does not extend up to very large parent 
polydispersities but terminates at a consolute point lo- 
cated dX a = 0.373 ± 0.001. This means that isotropic- 
nematic-nematic triphasic equilibria are only expected 
to occur in a small interval of parent polydispersities, 
namely 0.264 < cr < 0.373. However, it should be men- 
tioned that we have not been able to determine the cor- 
responding binodal points, which mark the actual on- 
set of three-phase separation from the isotropic-nematic 
biphasic equilibria. Considering the results for the ap- 
proximate Onsager 7'2-model, as numerically analyzed by 
Speranza and SoUich 18], we expect these triphasic equi- 
libria to be limited to a very small density interval across 
the coexistence region, which makes it very hard to ob- 
serve the phenomenon in experiment. Another problem 
is that the fraction of the system volume occupied by the 
nematic phases is predicted to be at the most 0.1 %, so 
that it will be very difficult to distinguish (or even detect) 
the two different nematics. Therefore we must conclude 
that, although the log-normal distribution contains suf- 
ficiently long rods to induce a demixing of the nematic 
phase, the (mole) fraction of these rods is too small to 
give rise to an observable fraction of the demixed nematic 
phases. 

Experimentally, phenomena such as a broadening of 
the biphasic region and a fractionation effect have been 
observed unequivocally in a number of experimental 
studies m 111 m llil. Observations of triphasic 
isotropic-nematic-nematic equilibria were however only 
reported for systems whose length distributions appear 
to be more or less bimodal rather than unimodal. These 
bimodal shapes were either accomplished deliberately by 
mixing species with different lengths, as done by Itou 
|26| with semiflexible schyzophyllan rods or caused by 
the presence of large aggregates as found by Kajiwara et 
al. \2i\ in systems of rigid imogolite rods. Buining et al. 
|25| observed the formation of second nematic phase in 
systems of polymer-coated (hard) boehmite rods, albeit 
a long time after the two-phase isotropic-nematic phase 
separation had finished. Also there, the triphasic demix- 
ing is probably due to the presence of a small number of 
very long rods or aggregates present in the system caus- 
ing the length distribution to be (slightly) bimodal [2^ . 

The experimental results therefore seem to indicate 
that having a three-phase separation with observable 
fractions of all phases requires some degree of bimodal- 
ity (with a sufficiently large length ratio between the 
short and long species) in the parental length distribu- 
tion. Hence, an intriguing issue left open for future in- 
vestigation is the question how the triphasic phenomenon 
predicted for the log-normal distributions change for a 
parent distribution with a slightly bimodal shape (e.g. 
with a small second peak just below Zmax)- In particular, 
a bimodal distribution may give rise to enhanced frac- 
tionation behavior and more pronounced triphasic equi- 
libria than predicted for the unimodal log-normal form. 



Appendix A: Numerical procedure 



The self-consistency equations for a{l) and /io^'(Z), 
Eqs. ^ and (|14|) . were solved using a numerical grid 
of lengths. The iterative scheme we used is analogous 
to the one described by Herzfeld et al. |29j for com- 
puting the numerically exact equilibrium ODF of the 
monodisperse Onsager model. The /—interval [Imim ^max] 
was discretized into N (not necessarily equal) parts. The 
mesh size must be chosen very carefully, particularly for 
large Zmax, because the distributions become consider- 
ably peaked at low polydispersities. Therefore, for par- 
ent polydispersities lower than approximately cr = 0.25 
we chose to divide the integration interval into three 
regimes. The interval [1 — 8ct^,1 -f- Scr^] in the vicin- 
ity of the peak, where the distribution changes rapidly, 
was discretized into equal parts, and the intervals 
[^min, 1 — 8cr^] and [1 + 8(7^, Zmax], whcrc the distribution 
is generally much smoother, were both discretized into 
equal parts. For parent polydispersities larger than 
a = 0.25 the entire interval was discretized into N equal 
parts. It has been proved sufficient to use N = 150 in or- 
der to obtain quantitatively reliable results. However, for 
the calulation of the full phase split (Sec. III-C) a smaller 
mesh size {N = 50) was used to limit the computational 
burden. Increasing the number of mesh points leads to 
only marginally different results in this case while the 
calculation time increased dramatically. For small poly- 
dispersities we took Eqs. (|10f) and H15|) as initial guesses 
and the successive iteration was performed until the fol- 
lowing convergence criteria were satisfied 



max \a{ln+i) - < 10 
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AW(Wi)-/ii^)(?„) 



< 10"^ (35) 



After each iteration step 10% of the new solution had 
to be added to 90 % of the previous one for the next 
iteration step to ensure the convergence of the method. 

The iteration algorithm we figured out to calculate the 
phase equilibria, in particular for the isotropic cloud and 
shadow phases and the full phase split, can be described 
as follows. First, the corresponding equilibrium forms 
for a{l) and /icx^''(Z) were calculated for a given set of 
starting concentrations. These results were then put into 
the self-consistency equations for the cloud and shadow 
concentrations to obtain new values. These equations 
were obtained by recasting the self-consistency conditions 
for the moments (e.g. Eq. ]27\i ) in an iterative form such 
that the concentrations could be calculated by a simple 
fixed-point iteration. 

Finally, for the new concentrations, corresponding 
forms for d{l) and p.ix'^ (l) were computed and substituted 
again into the self-consistency equations. This procedure 
was repeated until all concentrations had converged to 
within 10~^. To ensure convergence of this iteration loop 
a damping percentage of 80 % was used, which means 
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that only 20 % of the change was retained at each itera- 
tion step. 



Appendix B: High-cutofF scaling results 

In this appendix we focus on the properties of the 
nematic shadow phase in the cutoff-dependent regime 
(cr > (Jt) for systems with infinite cutoff lengths. For the 
exact Onsager model, Speranza and SoUich [l9| made 
a detailed analysis of these properties based upon the 
high-density (and high-cutoff) scaling forms of the exact 
nematic ODF. Here, we will not reproduce the analy- 
sis but merely show that the scaling form for the length 
distribution in the nematic shadow phase in the limit 
^max — > oo is analogous to the one obtained in Ref. [l9j |. 
Consequently, all scaling properties which follow from the 
Gaussian approximation must be exactly the same as the 
ones derived from the exact high-cutoff scaling results. 

The first step is to solve the coupled set of consistency 
equations, Eqs. JHl and (|14|l . In order to obtain analytic 
solutions for these nonlinear integral equations we exploit 
the fact that the nematic shadow phase is completely 
dominated by the longest rods in the system at a > at- 
When the cutoff length increases the length distribution 
in the nematic shadow will be more and more peaked at 
I = ^max- In the limit of infinite ^max it is therefore justi- 
fied to use the Ansatzp^^)(Z) — S{l — lma.x) which suggests 
an effectively monodisperse nematic shadow phase only 
containing the longest rods in the system. Substituting 
the delta-function in Eqs. lO and H14|l allows us to obtain 
asymptotic forms for the Gaussian variational parameter 
a{l) and the excess chemical potential jlix \l) of the 
nematic shadow phase. These expressions now read 



"(0 = ^max-^(V^ii 

I 



0, 



AixHO = — 23/Vl+^-H^Anax), 
'-max 



where is given by 




(36) 
(37) 

(38) 



Note that ^(//^max) always has a value between zero and 
unity. A close inspection of Eqs. I^Hjl and (p?7jl reveals 
that the variational parameter a{l) scales as a{l) oc If^^^^ 
whereas the excess chemical potential /i^(0 remains of 
the order 0(1) for all lengths. Using this in Eq. H26() 
we can write down a scaling expression for the length 
distribution in the nematic shadow, which in its general 
form reads 



pW(0=cst- 



CaxP*°^(Oexp 2cl°h~W{l/ln,,.) 



(39) 



(cf^)^ 

where W is a contribution of the order 0{1) 

m/lm..) = \nT{l/U.^) + /i(^'(VW). (40) 



Note that W attains its maximum W = 4 for I = Zmax- 
The scaling solution Eg. is completely analogous to 
the one found in Ref. [13 , the only differences being the 
exact form of W{l/lmax) and the constant est. However, 
since these contributions are both subleading in the limit 
^max oo they are irrelevant for the rest of the analy- 
sis and hence do not influence the scaling results. The 
similarity between the exact high density scaling analy- 
sis and the Gaussian approximation is also confirmed by 
the fact that both theories predict the same high-density 
scaling result for the nematic osmotic pressure, namely 
&/3H = 34^^ 

For the sake of completeness, let us now briefly outline 
the basic results of the large cutoff scaling analysis. For 
a comprehensive treatment of this subject the reader is 
referred to Ref. For a Schulz parent, we may use 

Eq. (1221 in Eq. $^ to obtain 

p(^) (0 = ICNl'tj' exp [el - >V(V?max)] , (41) 

where e = 2c^^ — (z + 1) and z ~ a^^ — 1. For very 
large I the exponent exp[£/] will be the dominating con- 
tribution. At a > at, the nematic shadow is supposed 
to be dominated by the longest rods and the distribu- 
tion p''^'(^) should therefore be an increasing function 
of length. This requires s to be positive and yields the 
condition Cq"'' > ^{z + I). Since the concentration of the 
cloud phase appears to decrease with increasing l^aK this 
then implies that the isotropic cloud curve (and hence the 
nematic shadow curve) has a finite lower bound for large 
cutoff lengths. These limiting solutions, for which e = 0, 
therefore read 



.(0) 



using Eq.(IlH|), with 



1 

" 2^' 
1 

6^ 

.(0) _ „(o) 
-1 — '-■0 ■ 



1 

2^ 



(42) 



These results are plotted 



in Fig. 1. To be consistent, let us now look for a solu- 
tion for the transition polydispersity at above which the 
nematic shadow phase for a Schulz distributed parent is 
completely dominated by the longest rods. We start with 
the concentration of the nematic shadow phase which is 
proportional to the integral over the normalized length 
distribution, i.e Cq^'' oc / p(^)(Od/. From Eq. gU we 
thus obtain 



oc L 



Z"exp[d- W(/Anax)]d/. (43) 



For the sake of convenience we have set l-adn equal to zero. 
Since the integrand is dominated by the exponent exp[£Z] 
for large I we may approximate the integral by bringing 
all slowly varying contributions in front of the integral 
sign and evaluating them at / = Zmax, which gives 



SN) 



cc l-tJFeM~Wil/l, 



oc I 



3 

max 



exp [ely, 



eL 



dl exp[el] 



(44) 
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The next step is to recast the latter equation into a scal- 
ing relation for e. Taking the logarithm on both sides of 
Eq. (O gives 



e oc (3 — z) 



In 



'max ^ ^^n 



o{i 



-1 ) 

max/ 



(45) 



From the known limits e^max ~^ oo and e J, for ^max — > 
oo we can deduce that e/max must increase more slowly 
than linearly with Z^ax- Consequently, the second and 
third terms in Eq. (|45l) are both subleading contributions 
so that we retain up to leading order 



e oc (3 — z) 



InL 



(46) 



which shows that e/max indeed increases logarithmically 
rather than linearly as we already anticipated. How- 
ever, in order to make this result fully self-consistent it 
is also required that z < 3 (and correspondingly a > 0.5) 
since e must be positive. This means that a — 0.5 is 
a lower bound for the cutoff-dependent regime in the 
limit /max ^ cxD. In other words, the transition from the 
low fractionation regime to the regime where the nematic 
shadow is completely dominated by the longest rods oc- 
curs exactly at at = 0.5 for Schulz parents with infinitely 
high cutoff lengths. 

We now turn to the average length in the nematic 
shadow phase which is related to the first moment density 
"^^^ cx Jlp^'^\l)dl. Analogously to Eq. gH) it follows 
^ exp[eZniax]/£/max and that the average 



that c[ OC /max 



length hence scales as (/) = c['^^ /cq^^ cx /max- Since the 
distribution in the nematic shadow phase is dominated 
by exp[£/] we expect that only rods whose lengths are 
of the order 0{l/e) smaller than /max contribute to the 
average length. We therefore can write 



(/) = /. 

and from Eq. 

(/) = /,, 



l-O 



1 - O 



eL 



In /,, 



(47) 



(48) 



This result shows that the average length in the nematic 
shadow in principle diverges for /max oo but the loga- 
rithmic correction causes the actual (/) to be significantly 
lower than /max- 

A similar treatment can be given for a log-normal par- 
ent distribution. However, the analysis for the log-normal 
case is even more involved and we will only present the 
basic results and refer to Ref. for details. First, 

the concentration of the isotropic cloud phase appears to 
have the following /max-dependence 



.(0) 



In^ L 



41n(H-(T2)/„ 



O 



In /n 



(49) 



which is crucially different from the Schulz case because 
the concentration of the cloud and shadow phases now 



tend to zero rather than approaching boundary values as 
in the Schulz case. Second, the average length can be 
shown to behave as 



(/) = /„ 



l-O 



1 



In^ L 



(50) 



Similarly to the Schulz distribution, the average length 
scales as (/) oc /max but the correction term is now consid- 
erably smaller which implies that the fractionation effect 
is much more pronounced for log-normal distributions at 
CT > Ut, as we already noticed by comparing Figs. 3 and 



Appendix C: Polydisperse spinodal instability 
criterion for the nematic phase 

In this appendix we show that the anomalous behavior 
of the coexistence pressure in Fig. (13) can be related 
to the local instability of the coexisting nematic phase. 
In particular, it is established that the (local) extrema 
in the coexistence pressure are directly related to the 
spinodal points where the nematic phase becomes locally 
(and hence globally) unstable. Let us first consider the 
Gibbs-Duhem equation, which in its general form is given 
by 



dP 



(51) 



with Pi = Ni/V. For a continuous density distribution 
c{l) this relation can be expressed in terms of functional 
derivatives 



6{bl3n) 
Sc{l') 



= / dlc{l) 



Sc{l') ■ 



(52) 



We now focus on the evolution of the length distribution 
of the nematic phase (in coexistence with an isotropic 
phase) across the coexistence region and denote this dis- 
tribution by c*(/). The shape of c*(/) depends uniquely 
on the nematic fraction 7, so that the curves in Figs. 12 
and 13 represent trajectories parametrized by 7. As to 
the osmotic pressure, let us now consider an infinitesimal 
change of the coexistence pressure (5(6/3H)* correspond- 
ing to an infinitesimally small dispacement Sc*{l) on the 
trajectory in Fig. (13). Using the Gibbs-Duhem relation 
Eq. (|52J), it then follows that 



s{bpuy = 



dl 



5c{l') 
die* {I) I dl 



Sc*{l') 

c'il') 

, 6f3f,{l) 



Sc{l') 



&*(/')- (53) 



;*(;') 



Since c* (/) > for all / the infinitesimal pressure change 
is only zero when 



dl' 



Sc{l)Sc{l') 



Sc*{l) = 0, 



(54) 



c*(0 
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which is precisely the spinodal instabihty criterion for 
polydisperse systems 30], albeit restricted to compo- 
sition fluctuations Sc*{l) which lie on the trajectories 
parametrized by 7. This result therefore shows that the 
local extrema in the coexisting pressure where 5(6/311)* = 
represent spinodal points indicating local instability of 



the coexisting nematic phase. The instability direction 
6c* (l) is defined as the change of the nematic length dis- 
tribution c(l) corresponding to an infinitesimal displace- 
ment ^7 on the trajectory in Fig. 13. The presence of the 
spinodal points for the nematic phase imply the existence 
of a isotropic-nematic-nematic triphasic region. 
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